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abstract 

Galaxy masses can be estimated by a variety of methods; each apphcable in different circumstances, 
and each suffering from different systematic uncertainties. Confirmation of results obtained by one 
technique with analysis by another is particularly important. Recent SDSS weak lensing measurements 
of the projected-mass correlation function reveal a linear relation between galaxy luminosities and the 
I depth of their dark matter halos (measured on 260 kpc scales). In this work we use an entirely 
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independent dynamical method to confirm these results. We begin by assembling a sample of 618 
^ ' relatively isolated host galaxies, surrounded by a total of 1225 substantially fainter satellites. We observe 
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the mean dynamical effect of these hosts on the motions of their satellites by assembling velocity difference 
histograms. Dividing the sample by host properties, we find significant variations in satellite velocity 
CO ' dispersion with host luminosity. We quantify these variations using a simple dynamical model, measuring 

. -^260" ^ dynamical mass within 260 h^^ kpc. The appropriateness of this mass reconstruction is checked 

by conducting a similar analysis within an N-body simulation. Comparison between the dynamical and 
lensing mass-to- light scalings shows reasonable agreement, providing some quantitative confirmation for 
' the lensing results. 

00 . Subject headings: galaxies: fundamental parameters, structure, mass function, halos 
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' 1. INTRODUCTION the lensing results requires an accurate understanding of 

fSJ , A/T 14- J ji-j-j-4-j-T,i 1 lens geometries. It is important to quantitatively confirm 

O*- ^ ' Mass clustered around galaxies distorts the local space- °, . , • , , ; i 

' ,. u J- 1- I, J. 1-1 J • these lensmg results by independent methods. 

— time, bending light rays which pass near them and m- t ;i ■ 1 1 • n • 1 i- 1 

--h j - r j-i-x 1 In this work we briefly examine galaxy mass-to- light 

• ducmg distortions m the images 01 more distant galax- ,. .,, . . , , , • , , 

r^ . • ru, 4- f 4-1, J- 4 4- 1,1 11- scalmgs with an independent dynamical method. We ex- 

les. Observation of these distortions enables weak lensing ° , r ; 1 1 • 1 i 1 

I ' 4 f 4-1, ■ 4- J 1 4-- f 4-- tract a sample of host galaxies surrounded by systems 

: measurements of the proiected-mass correlation function ^, , , , ; , , ,, m t-T- •; 1 m 

2 (PMCF: Brainerd, Blandford and Small (1996); Fischer f ^^^^^^^^f^ /'^"^ data collected by the Sloan Digital Sky 
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. shift Survey (Smith, Bernstein, lischer, & Jarvis 2001) ,i , • 1 ; ;i f ; , • , n 1 • i • 

k>( J 4-1, ci T^- -4 1 ci c 4. 1 method IS analogous to the stacking of galaxies used m 

■ and the Sloan Digital Sky Survey McKay, et al. 2002 : , . ^ r ^1 ti^ta^t^ ■ 1 1 1 4_ 1 

'.^ , , r, n^nr,\ 1 -ill- 1 .• lensing measurements of the PMCl" . It is closely related 

5h hereafter MOz reveal an approximately linear correlation , , , 1 r ■, 1 1 

rrt ■ 1 , 41, 1 • -4 f T T 1 • J 4U J 41, to the satellite galaxy analyses of Zaritsky, Smith, Irenk, 

. /T . between the luminosity of L > L* galaxies and the depth „ /-,„„„n t r, ■. , •; /iV^V^^n 

r- , , , , , , 4 4-1 11 • 1,- 1, 41, ^1 & White (1993) and Zaritsky & White (1994). 
of the dark matter potential wells m which they reside. ^ ' ■' ^ ' 



2. data, and the selection of host and satellite 

galaxies 



This close connection is not particularly surprising in a 
modern picture of galaxy formation (Seljak 2000; Guzik & 
Seljak 2001), though it requires rather uniform integrated 

star formation efficiency across halos of a range of masses. Data for this study are drawn from Sloan Digital Sky 

The weak lensing effects induced by galaxies are very Survey^^: a combined imaging and spectroscopic survey of 

small; peak distortions of background galaxies in the SDSS 10^ deg^ in the North Galactic Cap, and a smaller, deeper 

are only 0.5%. In addition, quantitative interpretation of region in the South. The imaging survey is carried out in 

drift-scan mode in five SDSS filters (u, g, r, i, z) to a lim- 
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iting magnitude of r<22.5 (Fukugita et al. 1996; Gunn et 
al. 1998). The spectroscopic survey targets a 'main' sam- 
ple of galaxies with r<17.8 and a median redshift of z^O.l 
(Strauss, et al. 2001). Velocity errors in the redshift sur- 
vey are ~30 km/s. For more details of the SDSS sec York 
et al. (2000) and Stoughton et al. (2002). For this study 
we extract from the SDSS main galaxy sample a catalog 
of 102,922 galaxies with accurately measured redshifts in 
a range from z=0.01 to z=0.2. 

Potential host galaxies must pass an isolation cut: they 
must be at least two times more luminous than any galaxy 
within a projected distance of 2 Mpc and a velocity 
window of ±1000 km/s. Around each such host, we se- 
lect as possible satellites all galaxies at least four times 
fainter than the host (1.5 magnitudes), within a projected 
distance of 500 kpc, and within a velocity window of 
±1000 km/s. Combining the requirement of faint satel- 
lites with the magnitude limit of the SDSS spectroscopic 
survey limits our host galaxies to r< 16.3. 

There are 13,737 galaxies with r<16.3 in the sample. 
Of these, 618 are both isolated and have at least one faint 
satellite. Most of the remainder fail the isolation cut. The 
total number of satellites is 1225. While the mean number 
of satellites around each host is only two, there are host 
galaxies with as many as 19 satellites. An object with this 
many satellites is clearly more a group with a dominant 
central galaxy than a galaxy with a few satellites. We re- 
move these by requiring that the combined luminosity of 
the satellites be less than the luminosity of the host. All 
but 5 of the host systems pass this cut. Absolute magni- 
tude distributions of host and satellite galaxies are shown 
in Figure 1 . Many still fainter satellites exist; we are pre- 
vented from identifying them by the magnitude limit of 
the SDSS spectroscopic survey. 

To reveal the mix of host galaxy types, we divide galax- 
ies into early (redder) and late (bluer) types using u-r color 
as described in Strateva, et al. (2001). The u-r distribu- 
tion for host and satellite galaxies is shown in the lower 
right panel of Figure 1. Overall the host galaxies are 75% 
early and 25% late types. It is important to stress that 
this division of early and late is luminosity dependent in 
a way which is strongly color dependent. In u, late types 
consitute 25% of galaxies at all luminosities, while in z, 
late types are > 50% of low luminosity galaxies, and < 
5% of high luminosity galaxies. This variation in the mor- 
phological mix with luminosity complicates comparison of 
these results to the M02 lensing results in the bluer bands, 
where early and late type galaxies have substantially dif- 
ferent mass-to-light ratios. 

3. DEPENDENCE OF SATELLITE DYNAMICS ON HOST 
LUMINOSITY 

To examine the relationship between satellite dynamics 
and host galaxy luminosity, we bin systems by host lumi- 
nosity, then examine the host-satellite velocity difference 
histograms for each luminosity bin. All luminosities are 
based on SDSS Petrosian magnitudes, k-corrected in the 
manner described in Blanton et al. (2001). An increase in 
the dispersion of satellite motions with host luminosity is 
clearly detected. 

We quantify this variation by fitting the velocity dif- 
ference histogram in each luminosity bin to the sum of a 



Gaussian plus a constant, where the constant term rep- 
resents the contribution due to random interlopers not 
physically associated with the host galaxy. The Gaussian 
widths fit to these distributions (ar ) are direct measures of 
the velocity dispersion characteristic of satellites for host 
galaxies in each luminosity bin. While models of satellite 
motions (White & Zaritsky 1992) suggest that the line of 
sight velocity dispersion of satellites may be non-Gaussian, 
our measured profiles are well fit by these models. Errors 
quoted for these Gaussian fits include the degeneracy be- 
tween peak height and width. 

4. MASS MODELING 

To compare these results to the weak lensing results of 
M02, we need to determine model masses for each lumi- 
nosity bin. The M02 lensing results quantify mass-to-light 
scalings by measuring M260; the mass inferred by fitting a 
singular isothermal sphere model to the observed density 
contrast within a radius of 260 kpc. We calculate here 
a very similar parameter, M2QQ . 

To derive this parameter we adopt a Jeans approach. 
This approach is unlikely to represent the situation in de- 
tail, as the dynamical times for satellites at projected sep- 
arations of 500 kpc approach the Hubble time. We 
address this concern briefly by reference to N-body simu- 
lations in the following section. We begin with the relation 
(Binney and Tremaine 1987): 
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where r is the radius at which we measure the mass, v"^ is 
the average radial velocity of satellites squared, v{r) is the 
number density of satellites as a function of radius, and (3 
is the velocity anisotropy 
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By measiiring the projected number density of satellites as 
a function of radius we determine v{r) cx r~^'^. Measure- 
ments of cTr in apertures from 100 to 500 kpc show 
no significant variation, implying that the second term in 
the parentheses of Equation 1 is zero. We assume that the 
velocity anisotropy [3~Q. If this is incorrect, the masses 
determined here will be biased by a factor of ~(1.0 - /?). 
Finally, we assume that our measured cr^, the mean square 
line of site velocity of all satellites within a projected dis- 
tance of 500 kpc of the host, is equal to v^. These 
are significant assumptions, which we test below by simu- 
lation. 

With these assumptions, our mass estimator reduces to 
the simple form: 
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where r is the 260 h ^ kpc radius to which we integrate the 
mass, and v"^ is the one dimensional RMS satellite velocity. 

5. COMPARISON TO GIF SIMULATIONS 

A variety of assumptions were made to arrive at Equa- 
tion 3. It is possible that one or more of these assumptions 
introduces systematic error. In addition, it is desirable 
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to relate the mass esimator M2QQ to theoretically favored 
quantities like M2oo- As a first test, we conduct a par- 
allel analysis within a simulated universe, extracting the 
observables used in our analysis from the simulation. By 
reconstructing masses in the simulation using MggJ^, and 
comparing these to -M200, it is possible to both test the 
validity of Mjgg^and understand its relationship to Af200- 

To conduct this test we turn to the GIF simulations 
of Kauffman, et al. (1999). GIF combines N-body simu- 
lations of the evolution of dark matter clustering with a 
semianalytic galaxy modeling scheme, including gas cool- 
ing, star formation, feedback from supcrnovae, and galaxy 
merging. The end result is a catalog of galaxies with B, V, 
R, I, and K magnitudes and stellar masses. Peculiar veloc- 
ities for each galaxy are determined from the N-body simu- 
lations. As a result, these velocities have the full structure 
we might expect to find in the data. The most important 
limitation of the GIF simulations for this comparison is 
their relatively low mass resolution. They resolve only the 
most massive and luminous host galaxies. 

To test om analysis of SDSS satellite motions, we con- 
duct a parallel analysis within the GIF simulations. We 
define as isolated all galaxies which are at least two times 
brighter than any other galaxy within their halo. We then 
search these for cases with satellites at least four times 
fainter. We use the velocity differences between hosts and 
satellites to populate velocity difference histograms, and 
sort these by host luminosity, in the same manner used 
in the SDSS satellite study. As a first check of Equa- 
tion 3, we use the simulations to check our assumptions 
about (3 and the relationship between and v^. We calcu- 
late the velocity anisotropy from the simulation, and find 
/3 = 0.06 ± 0.03. We also find the LOS velocity dispersion 
cr^ = 1.03 ± 0.03 X v-^:. While the GIF simulations are not 
a perfect match to our satellite sample, this provides some 
support for the assumptions we made in deriving Equation 
3. 

A more direct approach is to confirm that the data and 
simulations show a consistent relationship between the 
observables: ar and luminosity. This comparison of the 
SDSS and GIF ar vs. luminosity measurements is shown 
in Figure 2. The SDSS luminosity measurements are in 
i, and the GIF luminosities in Johnson I, but as both are 
converted to solar luminosities the comparison remains ap- 
propriate. Both the SDSS and GIF results are consistent 
with a simple model in which ar oc 

Since the GIF simulations provide values for the M20Q 
associated with each halo, we can compare measures 
of M2QQ derived from galaxy velocities to M200 within 
the simulations. This comparison suggests that M2go"~ 
0.7xM2oo- It is essential to note that this relation is de- 
termined only for the most luminous galaxies, and over 
only a factor of four in luminosity. While this compari- 
son is limited, it provides some confidence that Equation 
3 sensibly relates satellite velocities to masses. A more de- 
tailed understanding awaits simulations with higher mass 
resolution, in which more direct comparisons can be made. 

6. MASS-TO-LIGHT SCALINGS AND COMPARISONS TO 

LENSING RESULTS 

The measurements of M2QQ derived from satellite mo- 
tions are shown in Figure 3. The relationship between 



^260 light in each passband is well fit by a single 

power law. As observed in M02, these relations are con- 
sistent with a power law index of one (constant M/ L) in 
all bands except it, where a flatter relation, with a best fit 
power law index ^0.6 is observed. 

To compare these results to the M02 lensing results, we 
fix the power law index to one, and fit the M^qq vs. lumi- 
nosity data in each band to obtain values for M / L. Values 
for these M/ L ratios are given in Figure 3. The dynamical 
M/L values are consistent with the lensing values at the 
la level in most bands. They differ most strongly in u and 
g, where M02 suggests results will be very sensitive to the 
mix of host types. 

The masses of L* galaxy halos are proportional to light, 
and M/L values vary strongly with color. More detailed 
quantitative comparison of the lensing and dynamical re- 
sults will only be possible when we can measure more iden- 
tical samples of hosts and lenses. This is difficult. The 
magnitude limit of the SDSS spectroscopy biases dynam- 
ical studies toward hosts at z<0.05, while lens geometries 
prefer lenses at z^0.15. Still, the M/L values derived by 
independent methods agree reasonably, especially in the 
redder bands, where differences between the host and lens 
samples are probably less important. 

7. CONCLUSIONS 

Recent SDSS weak lensing results (M02) revealed a lin- 
ear relation between galaxy luminosity and mass on halo 
scales. To quantitatively test this result, we have measured 
the relationship between galaxy luminosity and satellite 
dynamics for a large sample of reasonably isolated host 
galaxies. We observe a highly significant increase in satel- 
lite velocity with host luminosity. 

To make direct comparisons to the weak lensing results, 
we apply a simple mass estimator, M!^^^ , which is closely 
analogous to the mass modeling used in M02. A first or- 
der test of the validity of this model is made by comput- 
ing it within the GIF simulations. The relationship seen 
here between Mjgg" and luminosity matches reasonably the 
relation seen by entirely independent weak lensing meth- 
ods. This confirms the essential conclusions reached in 
M02. The luminous and dark components of galaxies 
are strongly coupled on ^200 kpc scales. 

It is unknown whether the scaling relations observed 
here for super-L* galaxies apply at lower luminosity. Per- 
haps the most important future extension of these studies 
will be to less luminous galaxies. Additional comparison 
of lensing and dynamical mass estimates will require mea- 
surement of more closely related host and lens samples, 
and more detailed study of the extraction of mass from 
the observable PMCF and satellite velocity structure. 
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of The University of Chicago, Fermilab, the Institute for 
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tory, the Max-Planck- Institute for Astronomy (MPIA), 
the Max-Planck- Institute for Astrophysics (MPA), New 
Mexico State University, Princeton University, the United 
States Naval Observatory, and the University of Washing- 
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Fig. 1. — The four panels here provide an overview of the host and satelUte populations selected by this method. The upper left panel is 
a histogram of the number of satellites per host. The solid line is the full sample, the dashed (dotted) line is for those systems where the 
host galaxy is more (less) luminous than the sum of its satellites. The upper right panel shows the absolute r magnitudes of the hosts and 
satellites in bins of 0.1 mag (using Ho = lOOhkms^^ Mpc^^ throughout). The lower left panel shows the distribution of magnitude differences 
between hosts and satellites. Finally, the lower right shows the distribution of u-r colors for hosts and satellites. The vertical dashed line in 
this figure shows the approximate dividing line between early (redder) and late (bluer) type galaxies discussed in Strateva, et al. (2001). 
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Fig. 2. — This figure compares measurements of satellite velocity dispersion vs. luminosity for real SDSS satellites (diamonds) and simulated 
GIF satellites (boxes). The two measures are well represented by a single relationship ar oc L"-^. Vertical error bars represent the uncertainties 
in determination of the satellite velocity dispersion. Horizontal error bars represent the range of host luminosities for each bin. The results 
from GIF simulations allow us to test the applicability of our mass estimator. The similarity of the simulations and observations at this 
observable level give some confidence in the comparison. 
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Fig. 3. — This figure shows the relationship between MjgJ* and luminosity in each of the five SDSS bands. In each plot the data points 
are the MjgQ* estimates for each luminosity bin. In the g plot, an additional point has been added giving the approximate mass derived from 
satellite studies of Zaritsky & White (1994), scaled to a radius of 260 kpc assuming isothermality. The vertical arrow in each plot marks 
the luminosity of an L« galaxy in each band. The solid line in each plot represents the best fit lensing M/L for from M02. The dashed lines 
represent the best fit constant M/L model from these dynamical measurements. Note that all five figures are on the same scale. Most of the 
variation of mass with luminosity is seen for L > L* galaxies. 



